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ABSTRACT 

The broad-band spectral distributions of non-thermal sources, such as those of several known blazars, 
are well described by a log-parabolic fit. The second degree term in these fits measures the curvature 
in the spectrum. In this paper we investigate whether the curvature parameter observed in the spectra 
of the synchrotron emission can be used as a fingerprint of stochastic acceleration. 

As a first approach we use the multiplicative Central Limit theorem to show how fluctuations 
in the energy gain result in the broadening of the spectral shape, introducing a curvature into the 
energy distribution. Then, by means of a Monte-Carlo description, we investigate how the curvature 
produced in the electron distribution is linked to the diffusion in momentum space. To get a more 
generic description of the problem we turn to the diffusion equation in momentum space. We first 
study some "standard" scenarios, in order to understand the conditions that make the curvature 
in the spectra significant, and the relevance of cooling during the acceleration process. We try to 
quantify the correlation between the curvature and the diffusive process in the pre-equilibrium stage, 
and investigate how the transition between the Klein-Nishina and the Thompson regime, in Inverse 
Compton cooling, determine the curvature in the distribution at equilibrium. We apply these results 
to some observed trends, such as the anticorrelation between the peak energy and the curvature term 
observed in the spectra of Mrk 421, and a sample of BL Lac objects whose synchrotron emission peaks 
at X-ray energies 

Subject headings: acceleration mechanisms: statistical - galaxies: active - galaxies: BL Lacertae objects 
- galaxies: BL Lacertae objects (Mrk 421, Mrk 501, 1H 1426+428, 1ES 1959+650, 
Mrk 180, PKS 0548-322) 



1. INTRODUCTION 

A defining feature of the non-thermal emission from 
different types of galactic and extragalactic sources is 
that their spectra are described by a power-law (PL) 
over a broad photon energy range. In several sources, 
however, their spectra show significant curvature that is 
typically milder than that expect ed from an ex p onen- 
tial c ut-off. In previous papers iMassaro et al.l (|2004l 
120061 ) discussed the curvature observed in the broad 
band X-ray spectra of the two well known HBL (High- 
energy peaked BL Lac) objects Mkn 421 and Mkn 501. 
The basic idea was that this curvature was not sim- 
ply the result of radiative cooling of high energy elec- 
trons, responsible of the synchrotron and inverse Comp- 
ton emission, but that it was essentially relate d to the 
acceleration mechanism. IMassaro et al.l (|2004T ) showed 
that curved spectral distributions, in particular log- 
parabolic (i.e. log-normal) ones, develop when the ac- 
celeration probability is a decreasing function of the 
electron energy. In subsequent works, through the 
analysis of a large collection of X - ray observa t ions o f 
Mkn 421, ITramacere et al.1 (|2007| k iTramacerd (|2007f ): 

Electronic address: andrea.tramacere@unige.ch 



ITramacere et al.l (^009) pointed out that the observed 
anticorrelation between the peak energy and the curva- 
ture measured in the synchrotron Spectral Energy Dis- 
tribution (SED), could be used as a clear signature of a 
stochastic component in the acceleration process. Very 
recently, the log-parabolic law has been also applied to 
describe the spectra l distribution and evolutio n of some 
Gamma-ray bursts (jMassaro fc Grindlavll201ll ). 

The principal aim of the present paper is to investi- 
gate this scenario by extracting information on the ac- 
celeration processes using the curvature parameter mea- 
sured in the observed synchrotron and Inverse Compton 
(IC) spectra of Synchro- Self Compton (SSC) sources. We 
study the conditions in which the energy distributions of 
electrons, resulting from stochastic acceleration, can be 
approximated by a log-parabolic law and how its cur- 
vature evolves during their acceleration, and the role of 
IC cooling. We compare predictions from our theoret- 
ical descriptions with the curved spectra of some HBL 
objects. 

In Sec. [2] we give an intuitive picture to take into ac- 
count the effect of random fluctuations in the energy gain 
of particles and the role these play in determining the 
spectral curvature, as a consequence of the multiplicative 
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central limit theorem, and compare these results with the 
analytical solution of the diffusion equation, in the "hard- 
spheres" approximation. In Sec. [3] and Sec. 0] we give 
a more physical description of the problem, using hrst a 
Monte Carlo approach, and secondly by solving numeri- 
cally the momentum diffusion equation. We discuss the 
evolution of the curvature in the electron distribution 
as a result of momentum-diffusion before equilibrium is 
reached, and the role that synchrotron and IC cooling 
processes play on reaching the equilibrium. In Sec. [SJ we 
study the peak energy, fluxes, and curvature, trends in 
the SED of both the synchrotron and IC emission, look- 
ing for the fingerprints of the stochastic component. In 
Sec. |51 we show how our results can reproduce the spec- 
tral trends observed in a some HBLs, in particular we 
investigate the relation between the peak energy and the 
curvature, and between the peak energy and the peak 
flux. The good agreement between predictions and ob- 
served trends, confirms that the stochastic acceleration 
mechanism can play an important role in the physics of 
the blazars' jets and other SSC sources. 

2. THE LOG PARABOLA ORIGIN: ANALYTICAL 
APPROACH 

2.1. Statistical description 

In the statistical picture, the change in energy of the 
particles at each acceleration step n s is expressed as 
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where 7 is the Lorentz factor of the particle and e is the 
fractional energy gain. We here investigate the role of 
fluctuations of e, on the spectral shape of the accelerated 
particles. With this aim in mind, we express the energy 
gain fluctuations as 



e = e - 



X 



(2) 



where the random variable x nas a probability density 
function with zero mean value ((%) = 0) and variance 
at, and s represents the systematic energy gain, that 
we treat as a non-random variable and the probability 
density function of e is defined on the range e > 0. The 
particle energy at step n s can be expressed as: 
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where 70 is the initial energy of the particle. This 
equation clearly shows that the final energy distribution 
(71(7) = dN (7) / dj) will result from the product of the 
random variables £j. The determination of an analytic 
expression for the distribution resulting from the multi- 
plication of generic random variable is not an easy task 
(Glen et al. 2004). Using the simplifying assumption 
that the particles are always accelerated, namely the ac- 
celeration probability, P a , is set to unity and applying 
the multiplicative case of the central limit theorem (e.g. 
Cowan 1998) it is possible to show that the particle en- 
ergies will be distributed as a log-normal law: 
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where Nq, is the total number of particles, fi = (In 7), 
al = er 2 (ln 7). We can determine these two quantities 



by taking the logarithm of Eq. |3l 
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where we have ignored the covariance terms since we are 
assuming the energy gain at each acceleration step being 
independent on the one at the previous step. Remember- 
ing that (x) — 0, a x = a e , and ignoring the 4-th order 
term, we can write: 



/i = ln (70) + n s 
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This equation shows that the variance increases lin- 
early with the number of acceleration steps and it is pro- 
portional to a £ 2 . Substituting /j and a 1 into Eq. [H 
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Hereafter we will consider decimal logarithms (log = 
logio, c e = l/log 10 e w 2.3), to make easier a compar- 
ison of the curvature results form this paper with those 
presented in observational papers. Taking the logarithm 
of Eq. [51 and substituting the parameters from Eq. [5] we 
obtain: 
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where K includes all the constant factors. This is a log- 
parabolic law with the curvature (2nd degree in log 7) 
coefficient given by: 

(10) 
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The interesting physical insight of this equation is that 
the curvature of the particle energy distribution is in- 
versely proportional to the acceleration steps (n s ), and 
to the variance of the energy gain (a 2 ). In the case of 
P a < 1, the distribution at step n s will be given by the 
convolution of different log-normal distributions for each 
acceleration step, with the distribution at n s broader 
than that at n s — 1 and conta ining fewer particles, as 
already noted in iPeacockl (|1981[ ) . 

Similar results are obtained considering a constant en- 
ergy gain but a fluctuating number of acceleration steps. 
Assuming that after a time t the probability distribution 
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for the number of steps undergone by a particle is given 
by a Poisson law, it is possible to show that the energy 
distribution follows a log-parabola whose curvature term 
depends on the inverse of the mean number of steps mul- 
tiplied by the duration of the acceleration process. 

2.2. Diffusion equation approach 

The above statistical description provides an intuitive 
link between the curvature in the energy distribution of 
accelerated particles and the presence of a randomiza- 
tion process, such as the dispersion in the energy gain or 
in the number of acceleration steps. However, this ap- 
proach does not give a complete physical description of 
the processes responsible for the systematic and stochas- 
tic energy gain, ignoring other physical processes, such 
as the radiative cooling and injection rates, or the accel- 
eration energy dependence, necessary to give a complete 
description of the particles energy distribution evolution. 
A physical self-consistent description of stochastic accel- 
eration in a time-dependent fashion, can be achieved 
through a kinetic equation approach. Employing the 
quasi-linear approxi mation with the inclusion of momen - 
tum diffusion term (jRamatvl fl979t iBecker et all 12006( 1. 
the equation governing the temporal evolution of n(j) 
is: 
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where D p (p/,t) is the momentum diffusion coefficient, 
DA{l,t) — (2/j)D p ('j,t) is the average energy change 
term resulting from the momentum-diffusion process, 
and S(j,t) — —C(-f,t) + A(j,t) is an extra term de- 
scribing systematic energy loss (C) and/or gain (A) , and 
Q{"f,t) is the injection term. In the standard diffusive 
shock acceleration scenario, there are several possibilities 
for which one can expect that energy gain fluctuations 
will occur, due to the momentum diffusion term. In par- 
ticular, for the case of a turbulent magnetized medium, 
the advection of particles towards the shock due to pitch 
angle scattering may be accompanied by stochastic mo- 
mentum diffusion mechanism. In this scenario, parti- 
cles embedded in a magnetic field with both an ordered 
(Bo) and turbulent (SB) component, exchange energy 
with resonant plasma waves, and the related diffusion 
coefficient is determined by the spe ctrum of the p lasma 
waves. Following the approach of Be cker et al.l (2006) 
we describe the energy distribution W(k) in terms of the 
wave number k — 2ir/\ with a power-law : 

w{k) = sj^ = s_Bpf(ky^ 
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with q — 2 for the "hard-sphere" spectrum, q = 5/3 for 
the Kolmogorov spectrum, and q = 3/2 for the Kraich- 
nan spectrum, the total energy density in the fluctuations 
being 

U SB = / W(k)dk . (13) 

Jk a 

Under these assumptions the mome ntum-diffusion coef- 
ficient reads (jO'Sullivan et al.ll2009H : 

>-«(fHA^r— <»> 
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where (3 a = Va/c and Va is the Alfven waves velocity, 
p g = pc/qB is the Larmor radius, and X ma x is the maxi- 
mum wavelength of the Alfven waves spectrum. The ac- 
celeration time for particles with Lorentz factor 7, whose 
Larmor radii resonate with one particular magnetic field 
turbulence length-scale, is dictated by the momentum 
diffusion coefficient (D p ) as, 
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The spatial diffusion coefficient relates to the momentum 
diffusion coefficient through the relation, D x D p ss p 2 fi\ 
(Skilling 1975), hence the escape time of the particles 
from the acceleration region of size R, depends on the 
spatial diffusion coefficient through the relation, 
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The coefficients in Eq. Q21 and their related time scales, 
can be expressed as a power-law in terms of the Lorentz 
factor (7): 
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of a time. 



(17) 
and Aq have the dimension of the inverse 
Analytical solutions of the diffusion equa- 
tion for relativistic electrons are frequently discu ssed i n 
the literature since the early work by Kardashcv (1962), 
in particular for the case of the "hard-sphere" approx- 
imation. Neglecting the S and T esc terms in Eq. Q21 
and using a mono-energetic and instantaneous injection 
(71(7, 0) = NpS('y ~ 7o))i the solution of the diffusion 
equation is (|Melroselll969t iKardashcv 1962): 
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ie. a log-parabolic distribution, whose curvature term is: 
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This result is fully consistent with that found in the sta- 
tistical description, indeed Eq. [T5] and Eq. [5] have the 
same functional form in both the statistical and in the 
diffusion equation scenario, with t playing the role of n s , 
D p o the role of the variance of the energy gain (a 2 ), and 
A p o the role of log e. Hence we can write: 
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It is interesting to note, that in the case of the "hard- 
sphere" approximation, the curvature term is simply 
dictated by the ratio of the diffusive acceleration time 
(tr>) to the evolution time (t). 
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3. NUMERICAL APPROACH: MONTE-CARLO 
SIMULATION WITH MAGNETIC TURBULENCE 

In this section we demonstrate explicitly how the in- 
troduction of energy fluctuations leads to curved spec- 
tral distributions of particles. This is carried out using a 
Monte-Carlo (MC) approach. 

In our simulations, we considered 10 5 particles injected 
into the system with a cold mono-energetic distribution 
of Lorentz factors, with 70 = 1. To compare these results 
with the ones presented in Sec. [2] we remind the reader 
that, in the MC approach, the duration of the accelera- 
tion process t is the equivalent of the number of acceler- 
ation step (n s ) used in the statistical picture, and that 
the probability of the particle to be up-scattered or down- 
scattered in the MC realizations, can be expressed in the 
statistical approach as P(e > 1) and P(e < 1), respec- 
tively. The scattering probability of the particles is dic- 
tated by the intensity of resonant waves in the turbulent 
magnetic power spectrum. As a working hypothesis we 
assume that particles interact with a turbulent magnetic 
field whose power spectrum is expressed by Eq. Q21 In 
each scattering, the particles have probability (1 + /3a)/2 
of being up-scattered, and probability (1 — /3a)/2 of be- 
ing down-scattered. The energy dispersion of the parti- 
cle due to resonant scattering with Alfven waves will be 
(AE 2 ) ex (Ej3x) 2 t, where E = m e c 2 ^. Using the very 
good approximation for the var iance of the prod uct of n 
uncorrelated random variables ()Goodmanl ll962): 



a 2 {nx i )=Tl{x i ) 2 T i 
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and plugging Eq. [5] into the equation above, we get: 

2 

(AE 2 ) oc (Ep A ) 2 t oc 7o^ 2 (ne 4 ) = (lo^) 2 n s % (22) 

since E is the particle energy at time t (namely step n s ), 
we have E 2 = (m e c 2 7o£ ns ) 2 , from which follows: 
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In the following two sections (Sec. 13. H and. Sec l3.2[) we 
study the consequences of the structure in the magnetic 
turbulence, on the evolution of the particle spectra, fol- 
lowing their stochastic acceleration in the turbulent field. 

3.1. Hard- Sphere Turbulence 

Under the "hard-sphere" approximation (q — 2), the 
spatial diffusion coefficient does not depend on the par- 
ticle energy, since the exponent of Eq. [15] is q — 2 = 0. 
Only three independent parameters exist in this descrip- 
tion: the scattering time, the escape time and the veloc- 
ity of the scatterers. The spectra are purely determined 
by how many scatterings have been able to occur, the 
velocity of the scatterer, and what fraction of the in- 
jected particles have escaped out of the acceleration re- 
gion. The scattering time relates to the spatial diffusion 
coefficient by i sca t ~ D x /c. Similarly, the resulting accel- 
eration time relates to the spatial diffusion coefficient by 
iacc ~ D x /j3 2 A c k, t A //3%. Thus, for "hard sphere" tur- 
bulence, the scattering and acceleration time scales are 
independent of the particle energy (since there is equal 
energy density of scatterers which particles of all energies 
may resonantly scatter with). 



The left-hand panel in Fig. Q] shows the resulting in- 
stantaneous evolution of spectra for the "hard sphere" 
turbulence. The log-parabolic shape is maintained along 
the entire acceleration process, as shown by the solid lines 
representing the fit of the MC distributions by means of 
the law in Eq. [9] The evolution of the curvature param- 
eter, obtained from the r in the log-parabolic fit, and 
plotted in Fig. [2] with the red dashed line, clearly shows 
the trend due to the momentum diffusion, in agreement 
with the prediction from Eq. Q1J] (blue line in the plot) 
demonstrating the connection between D p q, ^, and /?a- 

3.2. Soft turbulence spectra 

To account for the effects of turbulent magnetic field 
spectra softer than the "hard-sphere" case, we also con- 
sider acceleration in Kolmogorov and Kraichnan type 
turbulence. We have therefore to include a fourth pa- 
rameter in the MC simulation, in addition to the three 
considered above: the turbulent field spectral slope q (see 
Eq. Q2D. 

The right-hand plot in Fig. Q] shows the evolution of 
spectra for the "Kolmogor ov" turbulence case. Simi- 
lar spectra were obta i ned b y lLemoine fc Pelletierl (|2003h . 
and lO'Sullivan etail (|2009f ). who integrated the trajecto- 
ries of charged particles in a turbulent magnetic field em- 
bedded in a fluid. Th e results are c o mpar ed to the quasi- 
linear theory results, iBecker et al.l (J2006) (solid lines in 
Fig. [IJ. We can identify two phases in the temporal 
evolution. In the first phase, the spectral energy distri- 
butions are more symmetric, and the curvature evolves as 
in the q = 2 case, while in the second phase they develop 
a low-energy power-law tail. Fig. [2] shows that, for the 
Kolmogorov (green line) and the Kraichnan (black line) 
turbulence, r is systematically larger compared to the 
"hard-sphere" case (red line), and that for t > 2 x t acci r 
approaches to an asymptotic value (r ~ 1.2 and r ~ 1.5, 
for q = 5/3 and q = 3/2 respectively) ruled by the expo- 
nential cut-off in the equilibrium distribution. 

4. NUMERICAL APPROACH: DIFFUSION 
EQUATION WITH STOCHASTIC COMPONENT 
AND LOSSES 

Both MC approach and statistical description are able 
to explain the link between the curvature in the energy 
distribution of accelerated particles and the presence of 
a stochastic energy gain term. In order to incorporate a 
more complete description, taking into account the com- 
petition between radiative losses and acceleration, and 
its influence on the curvature, we use the diffusion equa- 
tion approach, already outlined in Sec 12.21 by inserting 
into Eq. [12] a cooling ter m for the s ynchrotron and IC 
radiative losses. Following lModerski et al.l (|2005T > we can 
write: 



17 synch \ 



\1IC\ 



AaTC ;l 2 U B = C ~? 2 U B 



3m e c 2 



(24) 



7 2 / fKN('i-feo)£onph{eo)deo = Coj 2 F KN ('y) 
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where U B = B 2 /8ir, is the energy density of the mag- 
netic field, eo = hvo/m e c 2 is the IC seed photon en- 
ergy in units of m e c 2 , n p h(eo) is the number density 



Stochastic acceleration and spectral curvature in blazars 



— t= 


=0.6 t 

acc 


- t= 


=2.0 1 

acc 


— t= 


=6.0 t 

acc 



I I I I I I I I 

Hard-Spehere turbulence 
instantaneous spectrum 








iii 




| ' 


i i i i i 


- 







— t=0.6t 




Kolmogorov turbulence 






- t=2.0 1 

acc 




instantaneous spectrum 


- 






— t=6.0 t 

acc 






- 










"la 






- 


so 
£ -1 


J / Tl 


|_\ 
\ \ 
|\ 


\ \ 

Us 
\ 

L-A 

1 \ 


- 








\ 1 \ 






T 






- 






i 


^ v 
\ * ^ 


_ 






i 

\ 
\ 


r ^ 


\ 
\ 


- 




- \ . 


i . 


1,,,, 


- 



log (y) 



Fig. 1. — Plots showing the Monte Carlo results. For comparison, the results obtained using an analytic description given in lBecker et al.l 
(2006) (dashed line, right panel) and a log-parabolic function (dashed line, left panel), are shown. 
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Fig. 2. — The curvature parameters of the energy distribu- 
tion of accelerated electrons in the shown in Fig [l] In the case of 
q = 2 (red line), the trend is consistent with the "hard-spheres" 
prediction (blue line). In the case of Kolmogorov (green line) and 
Kraichnan (black line) turbulence, the trend predicts larger values 
compared to the "hard-spheres" prediction, and r approaches an 
asymptotic value dictated by the exponential cut-off in the equi- 
librium distribution. 



of IC seed photons with the corresponding photon en- 
ergy density U p h — m e c 2 J tQn p h(to)dtQ. The fu nction 
fxN r esults from the analytical integration of the iJonesI 
(1968) Compton kernel, fully taking into account Klein- 
Nis hina (KN) effects for an isotropic seed photon field 
(see lModerski et aT1l2005l appendix C), and Fkn(j) rep- 
resents its convolution with the seed photon field. We 
remark that Fkn pl a ys a crucial role in the cooling pro- 
cess, depending both on the IC regime (Thomson (TH) 
limit for 4760 << 1, KN limit for 4760 >> 1), and on 
eoriph(eo) oc B 2 /R 2 . 

Since analytical solutions are possible only for a lim- 
ited number of cases, to follow the complex dependence of 
the IC cooling term on n p h(to) in a self-consistent way, 
we must solve the diffusion equation numerically. For 
this purpose we further developed the numerical code 
(jTramacere et al.1 l2009t iTramacerel I2007 1 ) used to com- 
pute numerically the synchrotron and IC emission, and 
introduced it into the numerical solution of the diffusion 
equation. In the numerical calculations, we adopted the 



method proposed by Chang fc Cooper (1970) and use d 
the numerical recipe given by iPark fc Petrosiarj (|1996h . 
This is a finite difference scheme based on the centered 
difference of the diffusive term, employing weighted dif- 
ferences for the advective term. We use a 5000 points 
energy grid over the range 1.0 ^ 7 ^ 10 9 , and a time grid 
is finely tuned to have a temporal mesh several orders of 
magnitude smaller than typical cooling and acceleration 
time scales. The results from our code were compared, 
when possible, with known analytical solutions, and al- 
ways found good agreement. 

4.1. Physical set-up: the relations between D p , and tr> 
with j m ax and R 

We study the evolution of n("f) and of the curvature 
term in an homogeneous spherical geometry, with radius 
R and an entangled coherent magnetic field B and a tur- 
bulent component 5B, in the two cases of impulsive and 
continuous injection with a quasi mono-energetic source 
function Q{^i n j,t) normalised to have a fixed energy in- 
put rate: 



Linj — Tr 71 "-*?" 



r yin 3 rn e c 2 Q{^ in: j,t)d^i n: j (erg /a) (25) 



In our approach we don't distinguish the acceleration 
region from the radiative one, and during the acceleration 
process we take into account both synchrotron and IC 
cooling. According to Eq. [TJJ to determine the order of 
magnitude of D p we assume 1 >> SB/B ~ 0.1 — 0.01 
and require Alfven waves to be at least mildly rclativistic, 
with Pa — 0.1 — 0.5, and their maximum wavelength 
to be much smaller than the accelerator size (X ma x < 
R) . To study the effect of IC cooling on the evolution of 
71(7), we consider two different sizes of the acceleration 
region, a compact one (R = 5 x 10 13 cm), and a larger 
one (R = 1 X 10 15 cm). With this choice of accelerator 
size we set \ max ~ 10 12 cm. We stress that the choice 
of Xmax constrains the accelerative upper limit through 
P g < X max leading to j max < {X max qB)/m e c 2 , since 
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Fig. 4. — Left panels: evolution of the particle spectrum with impulsive injection and no escape for the case of R = 1 X 10 cm 
and q = 2. Upper panels represent the temporal evolution of 71(7), lower panels represent the temporal evolution of 7 3 n(7). Solid lines 
represent the case of SSC cooling. Red and blue solid lines, represent the final state for B = 1.0 G and B = 0.1 G, respectively. Green 
solid lines represent the temporal evolution, for B = 0.1 G, with step of 0.8 X £j> The dashed lines represent the final stage in the case 
of only synchrotron cooling. The vertical dot-dashed lines represent the equilibrium energy in the case on only synchrotron cooling. Right 
panels: Evolution of the curvature as function of t/to ■ Upper panel: curvature r evaluated at -y p , for the case of SSC cooling (solid red 
and blue lines) and for the case of only synchrotron cooling (dashed red and blue lines). The solid green line represent the prediction from 
Eq. 1191 Lower panel: the same as in the upper panel, for the curvature r% p evaluated at 73 p (empty and filled circles) compared to the 
case of r (solid lines). 



particles with larger p g (hence larger 7) can't resonate 
with shorter wavelengths. Taking into account a coherent 
magnetic field of the order of 0.1 G, and X max s* 10 12 cm 
we found that the purely-accelerative efficiency limits the 
particle energy to jmax ^ 10 75 . In the left panel of Fig. 
we plot to , given by Eq. [T71 as a function of X m ax , 
for the case of q = 2, SB/ B— 0.1, and /3 J 4=0.5. In this 
case the acceleration time is energy independent, and for 
Xmax ~ 10 12 cm it will be of the order of to = l/D p o « 
10 4 s. In the case of q ^ 2, the acceleration will have an 
energy dependence given by Eq. I17[ as shown in the right 
panel of Fig. [3] for the case of q = 3/2. In this section, 
we focus on the evolution of the curvature as a function 
of the momentum-diffusion term, and therefore use only 
the accelerative contributions coming from the diffusion 
terms (D p (7), Da{^)), neglecting the systematic extra 



term ^(7). All the parameters and their numerical values 
are given in Tab. 1. 

4.2. Impulsive injection 

In the left panels of Fig. 01 and Fig. [S]we plot the evo- 
lution of energy distribution n(j,t) (upper panels) and 
of 7 3 n(7,i) (lower panels) in the case of the impulsive 
injection without escape, for q = 2, and for two values of 
ft lx 10 15 cm (Fig. SJ) and 5 x 10 13 cm (Fig. EJ). We 
inject a quasi- monoencrgctic electron distribution with 
-finj w 10. The 7 3 n(7,i) representation is useful to com- 
pare the results concerning 71(7) presented in this section, 
with those regarding the synchrotron emission presented 
in Sec. [S] We denote by j p the peak energy of 71.(7) an d 
by r the curvature evaluated by means of a log-parabolic 
best fit over a one decade-broad interval centered at 7 p . 
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FlG. 5. — Left panels: The same as in Fig. [4] for the case of R = 5 x 10 13 cm. Right panels: Upper panel: the same as in Fig. [4] for 
the case of R = 5 x 10 13 cm. Lower panel: the evolution of the curvature r^p, for the case R = 5 x 10 13 cm. 
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case of R = 10 15 cm and B = 1.0 G, Bottom, case of R = 10 13 cm and 



73 P and r% p represent the peak of r y 3 n( , y) and its curva- 
ture, respectively. In the right panels of Fig. @J and[5j wc 
report on the corresponding temporal evolutions of the 
curvatures under the effect of both momentum-diffusion 
and cooling terms. The solid black line corresponds to 
t = 0.2 x tacci where t acc = t£> is the acceleration time 
due to momentum-diffusion. As the time increases, the 
diffusion term acts on the distribution by means of both 
Da and D p . The effect of the latter is to make the dis- 
tribution broader. 

One can distinguish three phases: in the first one the 
energy of particles increases and the curvature parame- 
ter decreases following a law r oc t" 1 in agreement with 
the statistical scenario of Sec. [2] and with the Eq. \T9\ 
independent of the magnetic field strength (B = 1,0 G 
and B = 0.1 G) and of the source size, because the accel- 
erative contribution dominates over the radiative losses; 
in the second phase, the radiation losses become relevant 
and the distribution approaches the equilibrium with an 



increase of the curvature; in the third phase, the balance 
between acceleration and radiation losses is established 
and the curvature reaches a stable value. 

The equilibrium distribution reached through stochas- 
tic acceleration, is described by a relativistic Maxwellian 
(jStawarz fc Petrosianl [20081 ). 



71(7) ex 7 exp 



-1 /7 



/(9,-y) 



(26) 



where /(<?, 7) is a function depending on the exponent 
of the diffusion coefficient and on the cooling process, 
and 7 eq is the Lorentz factor that satisfies the condition 



leq 

tcooiil) = t acc (-f) and is given by 

1 



leq = 



tn 



(27) 



0C0 (Ub+Fkn^)) ' 

with tacc equal to the fastest acceleration time scale 
among tA,to,tDA- In the case of Compton dominated 
cooling we have j eq oc - — tJj-t — , whilst in the case of 
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Fig. 7. — Left panes: evolution of the particle spectrum with impulsive injection and no escape for the case of R = 1 X 10 15 cm, B = 1.0 
and q = 3/2. Since tr) is energy dependent, on the x-axis we plot the ratio t/tr)('yinj), where tuiciinj^) ls the diffusive acceleration time 
evaluated at the injection energy -fi n j. Green solid lines represent the temporal evolution, for B = 0.1 G, with step of 2.4 X trj(-fo)- Right 
panels: Evolution of the curvature r (upper) and r^ v (lower). 
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panel: Evolution of the particle spectrum for continuous injection, R= lx 10 15 cm, B = 1.0 G, and q = 2. Lower panel: 
curvature r^ p . 



strong KN regime, or in general for synchrotron domi- 
nated cooling, we have -f eq ex - t — ^. Using a power- law 
form for the acceleration terms, and in the case of only 
synchrotron losses (or any cooling process that can be 
expressed as a power- law function of 7), it is possible to 
give an analytic expression of f(q , 7) (|Katarzvhski et all 
120061 IStawarz fc Petrosiad 12003 ). The expectation for 
synchrotron and IC/TH cooling process, and for q = 2 
is f(q, 7) = 3 — q = 1. The curvature resulting from a 
log-parabolic fit over a decade centered on j p is r ~ 2.5, 
and r^ p ss 6.0 in the case of 73^ 

We first discuss the case of R = 10 15 cm (Fig. [4]), 
with only synchrotron cooling (dashed lines, left panels). 
In terms of behaviour, we note that for the larger value 
of B (1.0 G, red lines, right panels), the r-t trend de- 
parts from the purely accelerative one (r oc t , green 
line right panels) early (relative to the B = 0.1 G case, 
blue lines in the right panels). This happens because 
the synchrotron equilibrium energy (vertical dot-dashed 
lines, left panels) is lower in the case of B = 1.0 G. 



For both values of B, the final values of r are close to 
the synchrotron equilibrium value of « 2.5. When IC 
cooling is also taken into account, the final values of the 
curvature in 71(7) are r w 2.5 and r « 0.6, for B = 0.1 
G and B = 1.0, G respectively. This difference is due 
to the different IC cooling regimes for the two cases. 
To show clearly the complexity of the transition from 
the TH to the KN regime, and its dependence on R 
and B, in Fig. [6] we plot the ratio jic/jsynch. (solid 
lines), and Uph/Us (dashed lines), as a function of 7 and 
normalized to unity, for the case of q = 2, for the fi- 
nal step of the temporal evolution. As long as the ratio 
Uph/Us is close to jic/jsynch., electrons cool in the full 
TH regime, and (7(7) = Cqj 2 (Ub + U p h)- On the con- 
trary, when the electrons radiate in the full KN regime 
lie lisynch. « Uph/Us- In this case, due to the inef- 
ficient KN cooling regime we have -jsynch. >> jic, and 
the cooling term is dominated by the synchrotron com- 
ponent: C{"f) « Co7 2 c7s. In the intermediate cases, it's 
difficult to estimate analytically the ratio jic/jsynch.- 
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TABLE 1 
Parameters values adopted in the numerical solutions of the diffusion equation for the cases studied in Sec 4 





impulsive inj. 


cont. inj. 


R 


(cm) 


5 x 10 13 ,1 x 10 15 


- 


- 


- 


B 


(G) 


0.1, 1.0 


- 


- 


- 


l^inj 


(erg/s) 


10 39 


- 


10 37 


- 


g 




2 


3/2 


2 


3/2 


to = 1/Apo 


(«) 


1 x 10 4 


1 x 10 3 


1 x 10 4 


1 x 10 3 


T- ■ 


to 


100 


- 


1 x 10 4 


- 


1 esc 


(fl/c) 


oo 


- 


2 


- 


Duration 


to 


1 x 10 5 


- 


- 


- 


"Jinj 




10.0 


- 


10.0 


- 



For B = 1.0 G, the SSC equilibrium is reached at 7 « 
3 x 10 4 and the SSC cooling occurs between the KN and 
TH regimes (see top-right panel in Fig. [6]), hence the 
value of / is different from unity, as predicted for the 
case of full IC/TH or synchrotron cooling. When B = 0.1 
G, the equilibrium energy is 7 « 10 7 and electrons are 
in extreme KN cooling (see top- left panel in Fig. 0, 
synchrotron losses are much higher than those due to 
IC scattering, and again r reaches the previous value of 
w 2.5. It is also interesting to note the difference in 
the trends of r-t and r% p -t. In the latter case, the trend 
departs form the purely accelerative regime earlier (see 
FiglU right-lower panel) since the electrons with energies 
close to 73 P are more energetic than those close to "f p , and 
thus have much shorter cooling times. 

The results for the compact region (fl = 5x 10 13 cm) 
are plotted in Fig. [5] Considering that the injected elec- 
tron luminosity is the same (see Table 1), we expect a 
different response from the IC cooling, due to the higher 
photon density n p h(eo)- The r evolution for the syn- 
chrotron cooling case is similar to the previous case, while 
for the SSC emission, both for the case of B — 1.0 G and 
B = 0.1 G, the final value of r is about 2.5. This is due 
to the larger photon density which moves the IC scatter- 
ing into the TH regime also for the case of B = 0.1 G 
(compare bottom-left to top- left panel in Fig. [6]), hence 
n(j) approaches the solution of Eq. [26] in the case of 

/ = I- 

In Fig. [7] we show the temporal evolution for the 
case of q = 3/2 (R = 1.0 x 10 15 cm, B=0.1 G). In this 
case, contrary to the q = 2 case, the acceleration time 
to is energy dependent, hence we study the evolution 
of r as a function of t/toilfinj), where toijinj) is the 
diffusive acceleration time evaluated at the injection 
energy Tin,-. The energy dependence of tn affects 
the evolution of r, and the shape of the equilibrium 
distribution, indeed, the r-t, and r% p -t trends show 
different behaviour compared to the case of q = 2. The 
equilibrium curvature is reached for t > 1 x tD^ifij), 
and the two equilibrium curvature values are r rs 1.2 
and T3 P « 3.0, roughly half of those found for the case 
of q = 2, and in agreement with the result from the 
MC. We note that, the curvature obtained by means 
of a log-parabolic fit of Eq. [26] for the case q = 3/2 
(namely / = 1.5), is r w 3.7. Hence, both the MC and 
the numerical solution of the diffusion equation give a 
result different from that predicted by the analytical 
solution in Eq. [26] 



4.3. Continuous Injection 

The case of continuous injection (see Fig. [8]) is more 
complex. The distribution developes a low-energy power- 
law tail, but a log-parabolic bending, driven by the dif- 
fusion, is still present at high energies, hence we evalu- 
ate the curvature only at 73^, (i.e. the representation 
useful to compare to the synchrotron SED emission). 
Spectral curvatures are generally milder than the im- 
pulsive injection. In the left panel of Fig. [8] we plot 
the r-t trend both for the case of impulsive (red lines) 
and continuous (blue lines) injection, the curvature in 
the continuous injection case are systematically lower 
in the pre-equilibrium phases, and in the acceleration- 
dominated stage the trend is again consistent with the 
"hard-sphere" approximation and statistical approaches. 
The slope of the electron distribution in the power-law 
tail is rj 1.06, in good agreement with the predicted one 
~ 1 + t min-accl (2t esc ) = 1.075, consistent with the re- 
sults of lKatarzvhski et"aTI (|2006[ ). 

5. EVOLUTION OF THE SPECTRAL 
PARAMETERS OF SYNCHROTRON AND IC 
EMISSION 

The most relevant parameters describing the SED of 
SSC sources provided by observations are the peak en- 
ergies and curvatures of the synchrotron and IC compo- 
nents. We denote these curvature parameters by b s and 
b c , respectively, and by E B , E c and S s , S c , we denote the 
corresponding SED peak energies, and flux values. We 
use v s and v c to indicate the corresponding SED peak 
frequencies. In the following, we describe the results of 
the relations between these parameters assuming that 
electrons are injected into the acceleration region with a 
quasi mono-energetic spectrum with 7^ ss 10, and us- 
ing an injection time of 10 4 s. We use the same working 
hypothesis for the momentum diffusion coefficient as in 
Sec. 14. R and add a systematic acceleration time for the 
first order process tA = 1.5 x 10 3 s, in order to produce 
E s values ranging between optical and hard X-ray ener- 
gies. We set the radius of the region ati? = 2xl0 15 cm 
and the same duration for the injection and acceleration 
processes, namely 10 4 s. We varied the other parameters 
of the model, B, q, and D p o to verify how they affect the 
relation between the observable ones. All the parameters 
and their variation ranges are summarised in Tab. 2. 

A phenomenological approach, based on t he 6- 
function approximation (jTramacere et al.l I2007L 120091 : 
iMassaro et all 120061 I2004D . is useful to address the ex- 
pected relation between the curvature parameters and 
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Fig. 9. — Left panel: Evolution of synchrotron (black dashed lines) and IC (red dashed lines) SEDs, for the case of to = 1.5 X 10 4 
s and q = 2 (top panel), and for the case of toi^inj) ~ 6.3 X 10 s and q = 3/2 (bottom panel). All the other parameters as reported 
in Tab. 2. The solid lines represent the SEDs averaged overt the full simulation period, and the blue dashed lines (top panel) represent 
the SEDs corresponding to the transition from TH to KN regime. Right panel: The temporal evolution of b s (black squares) and b c (red 
squares) as a function of t/tn Q , for the case of q = 2 (top panel), and q = 3/2 (bottom panels) The cyan line (top panel) represents the b s 
trend predicted for the synchrotron emission in case of <5— approximation. The dashed lines (top panel) represent the PL best fit of both 
b s (purple) and b c (blue) trends. 



TABLE 2 

Parameters values adopted in the numerical solutions of 

the diffusion equation for the cases studied in sec. 5 



which implies 



parameter 


range 


R 


(cm) 


2 X 10 15 


B 


(G) 


[0.01 - 1.0] 


J-'inj 


(erg/s) 


10 38 


Q 




[3/2 - 2] 


t A 


to 


1.8 x 10 3 


t Do = 1/Dp 


to 


[1.5 - 25] x 10 4 


T ■ 


to 


10 4 


1 esc 


(R/c) 


2.0 


Duration 


to 


10 4 


"Jinj 




10.0 



their connections with the peak energies and flux val- 
ues. According to th e stand ard synchrotron theory (e.g. 
iRvbicki fc Lightmanl ([1986)), in the 8~ function approx- 
imation, the synchrotron SED peak value and the corre- 
sponding peak energy can be expressed by the following 
relations: 



S s (E s )cxn(j 3p hi p B 2 6 
E.ocyLBS. 



(28) 



S s ex (E s ) c 



(29) 



where a — 1.5 applies for changes of 73^ leaving constant 
n(73 P ), a = 2 for variations of B only, and a = 4 when 
the main driver is 6. For a log-parabolic shaped ^1(7) we 
have: 

log(73 P ) = log(7 P ) + J: (30) 



and, using the relation b s ra r/5 (jMassaro et al.l 12004). 
or, more precisely, from the analysis presented in Sec. 
[4~2l b s ~ r 3p /5. It follows: 



log( J B s )«21og(7 P ) + — 
56 s 

The relation between b s and E s is: 

a 



log(E s /E ) 



(31) 



(32) 



with a = 3/5 = 0.6 

The spectral properties of the IC emission are more 
complex, depending on the transition fr om the TH to 
the KN regime (see iMassaro et al.l [2006, for a detailed 
discussion). In the former case, the curvature is close 
to that of the synchrotron emission, but systematically 
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smaller due to the energy redistribution by the scatter- 
ing process. In the transition to the KN regime, the 
energy of IC photons will approach 7m e c 2 , hence the IC 
spectral shape will reflect that of the high-energy tail of 
71(7), and the curvature b c will be closer to that of the 
electrons. Then, provided the IC scattering happens in 
the TH regime, the trends involving b c are expected to be 
similar to those of 6 S , but showing systematically b c < b s . 
As the KN regime is approached, b c changes differently 
from 6 S , converging towards r. 

5.1. Temporal evolution ofb s , and b c 

We compute the evolution of b s and b c , as a function 
of the time, for the case of £_d = 1-5 x 10 4 s, B = 0.1 
G, and q = 2, using temporal mesh of 2 s. We plot in 
the top-left panel of Fig. [9] the instantaneous SEDs at 
steps of 200 s: the solid lines represent the synchrotron 
and IC SEDs averaged over the full duration of the ac- 
celeration process (10 4 s). As the time is increased, the 
peak energy of both the synchrotron and IC SEDs moves 
towards higher energies with a broadening of the spectral 
distribution. The corresponding evolution of curvature 
parameters is reported in the top-right panel: b s has a 
trend similar to that of the electron distribution, with 
b s oc (t/tr)o)~ a , and a ~ 0.6 ( for comparison the cyan 
solid line represents the r^p/5 trend, as predicted by the 
S 5— approximation). The trend of 6 C , as expected, is 
more complex because of the transition from TH to KN 
regime. For t/t acc < 0.4, it follows the same trend of b s 
but with systematically lower values. For t/t acc > 0.4, 
when the TH-KN transition occurs, b c increases with 
time, approaching towards the electron curvature r value. 
This transition starts for values of E s « 5x 10 -3 keV 
(v s w 10 14 Hz), and E c w 0.05 GeV (v c w 10 22 Hz); and 
the corresponding SEDs are plotted by blue thick-dashed 
lines in the left panel of Fig. [9] 

In the bottom panels of Fig. HI we show the case of 
q = 3/2. The synchrotron curvature quickly approaches 
the equilibrium value of b s « 0.6, consistent with the 
equilibrium value r% p rs 3.0 discussed in Sec. 14.21 In this 
case we do not observe the TH/KN transition in the IC 
curvature, since the lower values of E s and -E c , keep the 
IC scattering mainly in the TH regime. 

5.2. E s -S s and E s -b s as a function of D p q and q 

The other parameter affecting the evolution of the 
spectral distributions is the diffusion coefficient D p q 
(see Eq. 15) which we assume to vary in the range 
[1.5 x 10 4 ,2.4 x 10 5 ] s _1 , studying how the main spec- 
tral parameters change. In the left top panel of Fig. [TO] 
we plot averaged SEDs for each different value of D p q. 
The top- right panel shows the trend of b c vs. D p q. As 
expected, for larger values of D p o, the curvature mea- 
sured at the peak energy is smaller. The trend is de- 
scribed by a PL with an exponent of about —0.6 for 
D p o < 2 x 10~ 5 s _1 , and with an exponent of about 
-0.25 for D p q > 2 x 10~ 5 s _1 . This break clearly shows 
the transition between the TH and KN regime (marked 
by a vertical dashed line) , indeed it happens for the same 
values of D p q corresponding to the TH/KN transition in 
both the D p o-b c trend, and the E c -b c plot (occurring at 
E c « 1 GeV, see bottom-right panel in Fig. [TOj)- The 
break in the D P Q-b s trend happens when electrons radiat- 
ing at E s enter the KN cooling region, hence, due to the 



lower cooling level (compared to the TH cooling regime, 
on the left side of the vertical dashed line) , the curvature 
decreases. 

Blue filled circles in the top-left panel represent the 
peak positions for both SED components. For the syn- 
chrotron component, according to Eq. I29[ the exponent 
a in the case of 72(73) =const, should be 1.5, while the 
results of the computations give a = 0.6. This difference 
is due to the fact that we inject in the mono-energetic 
initial distribution, always the same total power that cor- 
respond to the same number of particles. When the peak 
energy increases the distribution becomes broader, im- 
plying that the same total number of particles is spread 
over a larger energy interval and the number of particles 
contributing to the synchrotron peak emission decreases. 
Consequently, the S s -E s trend gets softer compared to 
the predicted value of 1.5. 

We verified quantitatively this effect by computing the 
trend n(73 P ) vs 7 2 p , and found a power-law relation with 
an exponent of about 0.98, in nice agreement with the 
difference between the exponent of 1.5 and that resulting 
in our simulations. In the bottom panels of Fig. [10] we 
plot b s vs E s (left), and b c vs E c (right). The S c -E c 
relation can be fitted by a power-law (orange line, top- 
left panel in Fig. [TU)) with the same exponent of the 
E s -S s relation, as long as the IC scattering, at E c and 
above, happens in TH regime. When the KN suppression 
becomes relevant (green line, top-left panel in Fig. ITUl) . 
the exponent is larger, and is close to unity. 

The synchrotron trend (bottom-left panel in Fig. [TO]) 
clearly shows the expected anti-correlation between the 
peak energy and the spectral curvature, which is well fit 
by the function given in Eq. 1321 with a = 0.68, not very 
different from 0.6, obtained for the (5-function approxi- 
mation of the synchrotron emission, and assuming that 
71(7) has a purely log-parabolic shape. A simple power- 
law fit of the same points returns an exponent —0.14. 

We also investigate the role of q on the spectral evo- 
lution, setting its variation range to [3/2 — 2], i.e. from 
the Kraichnan to the "hard-sphere" case. The relations 
between the spectral parameters are very similar to those 
found in the previous case with S s oc E® 6 , and S c oc E®- 9 . 
Also in this case, the synchrotron component follows the 
expectation with a lower curvature for harder turbulence 
spectra, and the IC trend shows the transition from TH 
to KN regime. The power-law best fit of S(E S ) vs E s 
gives a — 0.88, larger than that obtained for the case of 
D p q. In fact, for values of q lower than 2, correspond- 
ing to less turbulence and hence diffusion, the curvature 
gets higher values and the peak energy lower values, com- 
pared to the "hard-sphere" case. The power-law fit for 
b s vs E s returns an exponent of —0.16, practically coin- 
cident with the previous one, indicating that the average 
properties of these parameters are the same in both the 
q and D p q cases. 



5.3. E Kr -S 



s,c ^s,c 



id E SjC -b StC as a function of B 

The magnetic field B drives the radiative losses, which 
affect the evolution of the spectral parameters. In Sec. 
HI we showed that different cooling conditions, and the 
transition from TH to KN, can determine very different 
values of 7 eg , for the same acceleration conditions. As- 
suming that the acceleration time scale is independent 
of the magnetic field, Eq. [27] shows that j eq oc 1/-B 2 , 
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Fig. 10. — Upper left panels: synchrotron (red lines), and IC (red lines) average SEDs for each different value of tr> in the range 
reported in Tab. 2, with q = 2. Blue points represent the position of Eg q and Sg^c- The purple, orange, and green line represent the 
PL best fit of the Eg-Ss and Eq-Sc trends. Upper right panel: b s and b c , for each average SED in the right panel, as a function of D v q. 
Dashed lines represent the PL best fit of the b-D p o trend. Lower left panel: the b s -E s trend obtained by means of a log-parabolic best fit 
of the averaged SEDs plotted in the upper right panel Lower left panel: same as in the lower left panel, for the b c -E c 



implying that, as long as B is small enough to result in 
73p << 7eg, the evolution of 71(7) around the peak value 
is dominated by the acceleration terms, while, for values 
of B resulting in 73^ > j eq the evolution obtains a no- 
table contribution due to cooling. In the top-left panel 
of Fig. [TTJ we plot the averaged SEDs. According to Eq. 
l28l and [29l the synchrotron peak value should scale as 
S s oc (E s ) 2 . Indeed, for values of B < 0.2 G we obtain an 
exponent equal to 2.04, very close to the value found with 
the <5— approximation. For higher values of the magnetic 
field, E s is anti-correlated with S s . This behaviour repre- 
sents a cooling signature due to the decreasing value 7 e9 
for increasing B values, with 7 eg getting closer to 73 P . 
This is confirmed both by the shape of the synchrotron 
SEDs and by the b s -B plot (top-right panel in Fig. [TTJ. 
Indeed, S SEDs for B > 0.2 G, exhibit an exponential de- 
cay, meaning that the distributions have reached, or are 
close to reaching, the equilibrium energy. Consistently 
with the S shape evolution, the b s -B relation shows an 
almost stable value of b s for B < 0.2 G, and an increas- 
ing trend for B > 0.2 G. This change, in both the S s -E s 
and b s -B trends, is interesting and can provide a use- 
ful phenomenological tool to understand the evolution 
of non-thermal sources. Another interesting feature is 
shown in the S c -E c plot: for B < 0.2 G the IC peak en- 
ergy is practically constant, as expected in the KN limit 



from the kinematical limit relating the scattered pho- 
tons energy to that of the electrons: hvic ~ jrn e c 2 . In 
fact, photons at energies « E c are produced in the KN 
regime, and for B < 0.2 G the electron peak energy 7 3p 
is constant, so E c must also be constant. For B > 0.2 
G, 73 P decreases because of cooling, and, accordingly, E c 
also decreases. This is another interesting test that can 
provide a probe for B driven flares evolving to the KN 
regime. The E s -b s plot in the bottom-left panel of Fig 
[Til confirms the cooling signature discussed above, show- 
ing b s uncorrelated with E s as long as 73^ << j eq , and 
an increasing value of b s with E s almost stable, when 

73p > leq- 

6. SPECTRAL EVOLUTION OF HIGH ENERGY 
FLARES OF BRIGHT HBL OBJECTS 

The previous considerations on the spectral evolution 
of SSC sources, in which high energy electrons are accel- 
erated in a relatively short timescale by stochastic pro- 
cesses, can be successfully applied to describe the be- 
haviour of some bright HBLs objects. These sources are, 
in fact, characterised by having the synchrotron peak in 
the UV/X-ray range and the IC peak in 7 rays up to 
TeV energies. Several flares, observed simultaneously in 
both these ranges, exhibited SEDs very well described by 
a log-parabolic law, whose parameters, particularly their 
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Fig. 11. — Same as in Fig. 1101 for different values of B in the range reported in Tab. 2. 



curvature, are estimated with high accuracy. A similar 
analysis for low-energy peaked BL Lac objects is much 
more difficult because the peak of their synchrotron com- 
ponent is typically in the infrared range and the avail- 
able simultaneous m ultifrcqucncy d ata are ext r emely 
few. iTramacere eTail (|2007l 2009) and lTramaceTel (|2007f ) 
pointed out that the observed anticorrelation between E s 
and b s in the synchrotron SED of Mkr 421, can provide 
a clear signature of a stochastic component in the ac- 
celeration process. In the same analysis, these authors 
pre sented also an interest ing correlation between E a and 
S s . iMassaro et al.l (|2008D found that the E s -b s and E a - 



S s trends hold also for a larger sample of eleven HBLs, 
making stronger the hypothesis that a common acceler- 
ative mechanism may drive such physical processes for 
this class of active galactic nuclei . To give a theoret- 
ical framework to these phenomenological relations, we 
try to reproduce both the E s -b s and E s -S s relations de- 
rived from the data of the aforementioned papers. In 
the following we will conside r the data of Mrk 421 from 
ITramacere et al.l ((2007, 2009!), collected o yer a period of 
13 ye ars, and of six HBL objects from IMassaro et al.l 
(2008): Mrk 180, Mrk 501, PKS 0548-322, PKS 1959- 
650, 1H 1426+428, covering a period of about 11 years 
and includin g both quiescent and flaring states. The 
sources from IMassaro et all ()2008l ) were chosen because 
the data are good enough to safely constrain both cur- 
vature and E a values, and because the observed varia- 
tions of the sample luminosity are compatible with the 



assumption to be driven by changes of E s . 

Following the analysis presented in Sec. [31 we consider 
two scenarios in which these trends are driven by the 
momentum-diffusion term. In the first case, the momen- 
tum diffusion changes because of variations of D p q, due 
to changes of 5B/B or Pa, but the turbulence spectrum 
(q = 2) remains stable. In the second scenario the turbu- 
lence spectrum is variable with q ranging in [3/2, 2]. We 
use the same method described in Sect. [5] to compute 
the averaged SEDs for each value of D p (or q) ; computa- 
tions are performed for three values of the magnetic field 
B = 0.05, 0.1, and 0.2 G. All the model parameters are 
summarised in Tab. 3. 

The comparison with the data can be affected by an 
observational bias due to the limited energy range of de- 
tectors. In fact, when the peak energy is close to the lim- 
its the curvature is not well estimated because one can 
use only a portion of the parabola below or above the 
peak. Generally, curvatures lower than the actual ones 
are obtained. The energy range [0.5, 100.0] keV is the 
typical spectral window covered by X-ray and hard-X- 
ray detectors. In our analysis we used this fixed window 
to take into account this possible bias in the observed 
data when E s is variable. 

6.1. E s -b s relation 

The E s — b s trend, and in particular the anticorrela- 
tion between these two observables parameters, is the 
strongest signature of a stochastic component in the ac- 
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Fig. 12. — Left panel: the E s -b s trend observed for the six HBLs in our sample. The dashed green lines represent the trend reproduced 
by stochastic acceleration model, for the parameters reported in Tab. 3, and for the D trend, the different lines corresponding to three 
different values of B reported in Tab. 3. The purple lines represent the trend obtained by fitting the numerically computed SED over a 
fixed spectral window in the range 0.5 — 100 keV. Right panel: the same as in the left panel for the case of the q trend. 



TABLE 3 

Parameters' values adopted in the numerical solutions of 

the diffusion equation to reproduce the observed trends 

of the hbls reported in sec. 6 
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celeration. 

In Fig. [T2]we report the scatter plot in the E s -b s plane 
for the six considered sources. The left panel reports 
the results obtained by changing the value of D p0 : the 
green dashed lines describe the trend resulting from a 
log-parabolic fit of the synchrotron SED over a decade 
in energy centered on E s ; the purple lines represent the 
same trend obtained by fitting log-parabola in the fixed 
spectral window [0.5, 100.0] keV. Both these trends are 
compatible with the data and track the predicted anti- 
correlation between E s and b s . Purple data, however, 
give a better description, hinting that the "window" ef- 
fect could be a real bias. Each of the three lines was 
computed for a different value of the magnetic field. It is 
remarkable that the variation of a single parameter, D p q 
can describe the observed behaviour. The dispersion in 
the data is relevant, and can be related to the variation 
of B (as partially recovered by numerical computation), 
or by different values of the beaming factor, R, and Li n j, 
during different flares, and for different objects. 

The dot-dashed tick line represents the best fit of the 
observed data by means of Eq. 1321 and returns a value 
of a rs 0.6, as expected from theoretical predictions for 
the case of the £— approximation, and pure log-parabolic 
electron distribution. This fitted line is also compatible 



with the numerical trend shown by the purple lines. Note 
that the observed curvature values are in the range [0.1, 
0.5], corresponding to r 3p ~ [0.5, 3.0]. According to the 
results presented in Sec. 14. 2\ the expected equilibrium 
curvature in the synchrotron emission, in the full KN or 
TH regime, and for q = 2, should be of r% p rs 6.0, and 
of 7"3p ~ 5.0 in the intermediate regime. In the case of 
q = 3/2, the equilibrium curvature should be r% p ~ 3.0. 
This is perhaps an interesting hint that, both in the flar- 
ing and the quiescent states, for q = 2, the distribution 
is always far from equilibrium. In the case of q = 3/2, 
only for E s < 1.5 keV the curvature are compatible with 
the equilibrium (r^p ~ 3.0, corresponding to b s ~ 0.6). 
For larger values of E s , we find again curvature well be- 
low the equilibrium value. These results provide a good 
constraint on the values of the magnetic field B < 0.1 G. 
The g-driven trend (right panel) is also compatible 
with the data, but for values of E 8 < 1 keV, the D p q- 
driven case seems to describe better the observed be- 
haviour, but any firm conclusion is not possible because 
of the dispersion of the data. 

6.2. E s -L s trend 

As a last benchmark for the stochastic acceleration 
model, we reproduce the observed correlation between 
E s and S s , which follows naturally from the variations 
of D p o and q. Considering that the redshifts of the six 
considered HBL objects are different, we prefer to use 
their peak luminosity L s = S s 4-7T D 2 L , where Dl is the 
luminosity distance . To account for the different jet 
power of sources, we considered two data subsets, and 
we assumed Li n j = 5x 10 39 erg/s for the first subset 
(top panels of Fig. [T3|) . and Li n j = 5 x 10 38 for second 
(bottom panels of Fig. [T3|) . In the left panels of Fig. Q21 
we report the D p $ driven trend, and in the right panels 
the q driven trend. Solid lines represent the trend ob- 
tained by deriving L s from the log-parabolic best fit of 
the numerically computed SEDs, centered on E s ; dashed 
lines are the trends obtained by fitting the numerical re- 
sults in the fixed energy window [0.5, 100] keV. 



We used a flat cosmology model with: Hq- 
km/s/MpC, £Vattcr=0.27, Sivacuum = 0.73 
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Fig. 13. — Left panels: the E s -L a trend observed for six HBLs in our sample, top panel corresponds to the case of Li n j = 5 X 10 39 erg/s, 
bottom panel corresponds to the case of Li n j = 5 X 10 38 . The solid black lines represent the trend reproduced by stochastic acceleration 
model, for the parameters reported in Tab. 3, and for the D trend, the different lines corresponding to three different values of B reported 
in Tab. 3. The dashed lines represent the trend obtained by fitting the numerically computed SED over a fixed spectral window in the 
range 0.5 — 100 keV. Right panels: the same as in the left panel for the case of the q trend. 



Both the results give a good description of the observed 
data, and their shapes are similar. Solid lines follow well 
a power-law with an exponent of about 0.6, while the 
windowed trends (dashed lines) show a break around 1 
keV and the exponent below this energy turns to about 
1.5. A similar break at the same energy, can be noticed 
in the points of M rk 421 in the E s -S s plot presented by 
iTramacere et a l. (2009), who found an exponent of ~ 1.1 
and of ~0.4 below and above 1 keV, respectively. This 
could again be an indication that the observed values are 
actually affected by the bias. 

7. DISCUSSION 

Broad band observations of non-thermal sources have 
shown that the spectral curvature at the peaks of their 
SEDs can now be measured with good accuracy. In this 
paper, we have presented, using different approaches, 
the relevance of these data for the understanding of the 
competition between statistical acceleration and radia- 
tion losses. First, using a simple statistical approach and 
Monte Carlo calculations, we have shown that the log- 
parabolic energy distribution of the relativistic electron 
is a good picture in the first phases before equilibrium is 
reached. In this case the curvature decreases with time 
and, therefore, for increasing peak energies. This evolu- 
tion is confirmed by numerical solutions of the diffusion 
equation taking properly into account both stochastich 
acceleration and radiative SSC cooling. The major re- 
sults can be summarised as follows. 

The evolution of the electron energy distributions (Sec. 
Q} shows that: 

• in the case of synchrotron and SSC cooling, and for 
all the values of B and R, as long as the distribution 
is far from equilibrium, the trend on r is dictated 
by Dp, and is well described by Eq. [19] 



• when the distributions approach equilibrium, the 
value of r is determined by the shape of the 
equilibrium distribution, which is a relativistic 
Maxwellian, with the sharpness of the cuf-off de- 
termined by both q and the IC cooling regime; 

• in the case of q — 2, and for equilibrium energies 
implying that IC cooling happens either in the TH 
regime or in the extreme KN regime (IC cooling 
negligible compared to the synchrotron one), the 
numerical solution of the diffusion equation follows 
the analytical prediction (/ = 1, that holds for any 
7 oc 7 2 ), and the corresponding equilibrium curva- 
ture is r 3p ?s 6.0 (b s ~ 1.2). In the case of q = 3/2 
the equilibrium curvature is r?, p f» 3.0 (b s ~ 0.6). 
These limiting values could be a useful observa- 
tional test to find cooling dominated flares with 
the distribution approaching to the equilibrium; 

• when cooling is in the intermediate regime between 
TH and KN and for the q = 2 case, the condi- 
tion / = 1 fails, and the end values of r decrease, 
strongly depending on the balance between Ub and 
the seed IC photon energy (JJ p h)\ numerical com- 
putations are necessary to evaluate the right value 
of r at equilibrium. 

The analysis of the spectral evolution of SSC emission 
(SecE]) shows that: 

• changes of D p q (or q) imply that the curvature and 
peak energy of the synchrotron emission are anti- 
correlated; the E s -b s trend can be phenomenologi- 
cally described by Eq. [32j 

• The E c -b c trend presents a clear signature of the 
transition from the TH to the KN regime. In par- 
ticular when the IC scattering approaches the KN 
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regime we observe a sharp change in the 6 C , with 
a positive correlation with E c , whilst in the TH 
regime the correlation is negative as in the case of 
the E s -b c ] 

• the magnetic field plays a relevant role on the cool- 
ing process, and B driven variations present rele- 
vant differences compared to those due to D p o (and 

In particular, for what concerns the B driven case, we 
note first that the E s -S s correlation follows the predic- 
tion of the synchrotron theory and shows the power-law 
relationship with E s ex (S , s )~ 20 . On the contrary, in the 
case of D p o and q changes, we find E s ex (S s )°- e . Another 
relevant difference in the B driven case is the evolution 
of S c . For the case of D p q and q driven trends S c re- 
lates to E c through a power-law with exponent of about 
[0.7-0.8]. On the contrary, for the B driven case with 
IC scattering in the full KN regime, the value of E c is 
almost constant and uncorrelated with S c (see Fig. ITT]) . 
due to the kinematic limit of the KN regime. E c starts 
to decrease when B is enough large to make dominant 
the cooling process. This is an interesting signature that 
could be easily checked in the observed data. 

The comparison of the E s -b s and E s -S s trends, ob- 



tained through several X-ray observations of six HBL 
objects spanning a period of many years, with those pre- 
dicted by the stochastic acceleration model, shows very 
good agreement. We are able to reproduce these long- 
term behaviours, by changing the value of only one pa- 
rameter (D P Q or q). Interestingly, the E s -S s relation fol- 
lows naturally from that between E s and b s . This result 
is quite robust and hints at a common accelerative sce- 
nario acting in the jets of HBLs. 

As a last remar k , we note that very recently 
iMassaro fc Grindlavl (|2011h find also in the case of 
GRBs a E s -b s trend, similar to that observed in the case 
of HBL objects. They measured values of the curvature 
up to 1.0, typically higher than in HBLs. It's interesting 
to note that the value of 1.0 is close to the limit of ~ 1.2, 
that we predict in the case of distributions approach- 
ing the equilibrium in either TH or KN regime, for q = 2. 
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